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Engineering of the cross-section shape and size of ultra-scaled Si nanowires (SiNWs) 
provides an attractive way for tuning their structural properties. The acoustic and 
optical phonon shifts of the free-standing circular, hexagonal, square and triangular 
SiNWs are calculated using a Modified Valence Force Field (MVFF) model. The 
acoustic phonon blue shift (acoustic hardening) and the optical phonon red shift 
(optical softening) show a strong dependence on the cross-section shape and size of 
the SiNWs. The triangular SiNWs have the least structural symmetry as revealed by 
the splitting of the degenerate flexural phonon modes and The show the minimum 
acoustic hardening and the maximum optical hardening. The acoustic hardening, 
in all SiNWs, is attributed to the decreasing difference in the vibrational energy 
distribution between the inner and the surface atoms with decreasing cross-section 
size. The optical softening is attributed to the reduced phonon group velocity and the 
localization of the vibrational energy density on the inner atoms. While the acoustic 
phonon shift shows a strong wire orientation dependence, the optical phonon softening 
is independent of wire orientation. 
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I. INTRODUCTION 



Silicon nanowires (SiNWs) have the potential to become the next generation solution for 
improved (i) Complementary-Metal-Oxide-Semiconductor (CMOS) transistors with higher 
channel drive currents^, (ii) thermoelectric (TE) devices^, (iii) rechargeable batteries^, (iv) 
compact non-volatile memories^, (v) solar cellsp and (vi) explosive and bio-sensors^. Ex- 
tensive application of SiNWs to various technologies requires to understand the structuraP, 
electricaP, thermaP^and opticaP properties which arise from the strong geometrical confine- 
ment, atomic positions, and increased surface-to-volume ratio (SVR). Especially the phonon 
spectra of SiNWs can provide a lot of vital information about their structural properties 
which can be optimized for maximum technological benefits in SiNWs. 

Engineering the Si nanowire structural properties based on their cross-section size has 
been studied extensively^^. However, the impact of cross-section shape has not been 
completely addressed yet. Experimentally, the impact of semiconductor (GaAs, GaP, ZnS, 
etc) nanowire shape on their structural properties has been observecP^'. Motivated by 
these experiments, we theoretically explore the effect of the cross-section geometry of ultra- 
scaled SiNWs on their structural properties like phonon shifts. We furthermore provide 
analytical expressions for the variation of these physical quantities with size and shape thus, 
enabling (i) physical insights, and (ii) a compact modeling representation of the structural 
properties. 

The determination of size and stability of SiNWs is crucial for their usage in different 
applications^ 2112 ^'. Raman spectroscopy provides a valuable, non-destructive method to 
analyze the physical properties of nanostructures (NS) like cross-section siz d 11 ' 121 , degree 
of crystallinity^, dimensionality of confinement! 11 * 21 *, etc. Raman spectroscopy utilizes the 
amount of phonon spectrum shift and line width broadening to obtain useful structural 
information^. 

The first satisfactory theoretical model for explaining Raman shifts and line width broad- 
ening was proposed by Richter et aP^and improved by Campbell and Fauchet 24 . These phe- 
nomenological continuum models are based on the choice of the phonon weighting function 
as well as the boundary conditions of the phonon vibration d 23124 1 which are chosen rather 
arbitrarily and hence do not provide correct phonon shift in nanostructures^. Arora et al. 
experimentally^ reported that phonon shifts are very different from bulk in nanostructures 
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with size less than 20nm. Tanaka et alP also pointed out the failure of continuum models 
to predict the Raman shift in nanostructures with size < 7.5nm due to the arbitrary param- 
eter sets used in continuum models. While continuum models can match the experimental 
Raman spectroscopy data using fitting parameters, they lack the predictive capability for 
estimating the phonon shifts and line widths in nanostructures with different shapes and 



compositions^"^. 

The evident inability of the continuum models to properly explain and predict the struc- 
tural properties in nanostructures^HZB demands for methods which can automatically take 
into account the effects of structural miniaturization like geometrical confinement, orienta- 
tion effects, cross-sectional shape and surface-to-volume ratio effects. 

Many efforts to explain the phonon shifts in nanostructures have been reported. They are 
based on models like the VFF approach! 13 * 28 *, partial density approach 25 , Bond Charge Model 
(BCM)Q3ES bond order modeP^, etc. In all these approaches the importance of taking into 
account the atomic positions and the atomic surrounding to understand the Raman shift 
has been stressed which most of the continuum models lac k 23 1 24 1 . 

We utilize a Modified Valence Force Field (MVFF) phonon modeP^^ to study the phys- 
ical properties of SiNWs. This MVFF model (i) is relatively straight forward to understand 
and implement, (ii) explains the experimental phonon and elastic properties of zinc-blende 
materials very we l | 28 l 3 ° l 2 1l33] ; anc i (jjj) j s computationally very efficient which allows for ex- 
panding the simulation domain easily to few million atoms^. However, frozen phonon models 
like the MVFF one do not integrate thermodynamics which forbids the correct calculation 
of many physical properties such as low temperature volume thermal expansion. Yet these 
models give a good insight into the physical processes governed by the lattice under ambient 
conditions. 

The paper is organized as follows. Section |TT] provides a brief introduction about the 
MVFF model for the calculation of phonons in SiNWs and about the method to extract the 



phonon shifts in these nanowires (Sec. II B). The results section (Sec. Ill) discusses the effect 



of cross-section shape and size on the phonon spectra (Sec. Ill A) and presents experimental 



benchmarking of the phonon shifts using the MVFF model (Sec. Ill B ). Furthermore, the 



trends of the acoustic phonon shift (Sec. Ill C ) and the optical phonon shift (Sec. HID) are 



discussed along with wire orientation dependence (Sec, HIE). Conclusions are provided in 
Section ED 
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II. THEORY AND APPROACH 



A. Phonon model 



In the MVFF model, the frequencies of the phonon modes are calculated from the forces 
acting on atoms induced by finite displacements of the atoms from their equilibrium positions 
in a crystaP^. To calculate the restoring force(F), first the total potential energy of the solid 
(U) is estimated. For the MVFF model, U is approximated as^, 



U 
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where Na, nn{i) and COPi represent the total number of atoms in one unit cell, the 
number of nearest neighbors for atom 'i', and the coplanar atom groups for atom 'i', re- 
spectively. The first two terms U^ s and represent the elastic energy coming from bond 
stretching and bending between the atoms connected to each other as given in the origi- 
nal Keating VFF modeP^. The terms U^l_ bs , Ul l s k _ bb and U^ bb represent the cross bond 
stretching®', cross bond bending-stretching®', and coplanar bond bending^ 1 interactions, re- 
spectively. The detailed procedure for obtaining the phonon spectra in Si bulk and NWs 
are provided in RefJ 31 * 32 l In the following sections the method to extract the acoustic and 
optical phonon shifts and the types of SiNWs used in this study are discussed. 



B. Phonon shift in Si nanowires 

The vibrational spectra of the nanostructures are strongly influenced by (i) the finite size 
of the nanostructures, (ii) the mismatch with the matrix material ( the surrounding atoms) or 
(iii) defects ^ * 20 * 22 ! The optical and the acoustic phonon energies are shifted compared to the 
bulk values (around q oc 0) in these geometrically confined structures. Raman spectroscopy 
reveals the decrease of the optical frequency with down-scaling of the nanostructures (called 
the optical softening or optical red shift )P332 anc j the increase of the acoustic frequency 
(called the acoustic hardening or acoustic blue shift )li 2 -12ZI. These phonon shifts are estimated 
using the following relation ^ 4 * 27 * 29 !, 
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A, ac/opt _ ac/opt ac/opt /n\ 

where is the bulk acoustic (optical) phonon frequency at the T point. The phonon 

shifts show the following functional dependence on the nanostructure cross-section size 

Ay-|11I27I2HI 

£*tf* = A ac/opt Q)^ pt [cm -i ]} (3) 

where a is the lattice constant of the material. For A>0 the phonon shift is positive (or 
blue shift) and for A<0 the shift is negative (or red shift )^. The exponent value 'd' indicates 
the dimensionality of the geometrical confinement. For ID nanowires/nanorods a variety of 
'd' values have been reported in the literature like 0.67^, 0.9^P, 1.44^, etc. 

Using a series of Raman spectroscopy measurement for a sample of nanostructures the 
size and the dimensionality of the nanostructures can be predicted 15 * 22 * 27 ! The phonon shifts 
observed in Raman spectroscopy can be predicted using the phonon spectrum. Atomistic 
phonon modeling correctly takes surface boundary conditions into account depending on the 
cross-section shape of the wire which reflects in the overall phonon spectrum. The surface 
and inner atomic movements vary with different cross-section shapes and therefore affects 
the phonon shifts. Equation [3] captures the nanostructure size effect in the exponent value 
'd' and the shape effect in the pre-factor 'A'. Thus, a combination of A and d can be useful 
for the determination of size and shape of the SiNWs. 

C. Si Nanowire details 

Four cross-section shapes for [100] oriented SiNWs have been considered in this study, (a) 
circular, (b) hexagonal, (c) square and (d) triangular cross-section shapes (Fig. [T]). SiNWs 
with [110] and [111] channel orientation studied here are of square cross-section only. The 
feature size is determined by the width parameter W. The value of W is varied from 2 to 6 
nm. The surface atoms are allowed to vibrate freely without any passivating species. The 
wires are still assumed to have a tetrahedral geometry. It has been shown that wires with 
diameter below 2nm tend to lose the tetrahedral structur e! 38 1 39 1 due to surface pressure and 
internal strain. 
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III. RESULTS AND DISCUSSION 



In this section the effect of cross-section on the phonon dispersion and phonon shifts 
of [100] SiNWs are discussed. There is a lack of proper channel orientation information in 
the literature. However Adu et alW^ reported [100] channel orientation due to which we 
chose [100] SiNWs for shape study. Furthermore, we also investigate the impact of channel 
orientation on the phonon shifts in square SiNWs. 

A. Phonon dispersion 

The phonon dispersion of SiNWs for all the cross-section shapes with [100] channel are 
shown in Fig. [2j While the full phonon dispersion extends up to ~65meV, phonon modes 
till 9 meV energy range are shown for clarity. All the wires exhibit two flexural branches^ 
(u q ^o oc q 2 ). These modes are double degenerate in all the shapes except the triangular 
wire. In the triangular wires these modes split due to the reduced structural symmetry 
(Fig. [2] d). The next two phonon branches in all the structures show co q ^o oc q. These 
branches determine the sound velocity in these structures^ (Table W). The triangular wire 
has the lowest longitudinal and transverse sound velocity (V sn d,i) whereas other shapes have 
very similar sound velocities. The splitting of the degenerate flexural modes and the reduced 
sound velocity show that triangular wires have the least structural symmetry. 

B. Experimental Benchmark 

Before using the MVFF model to calculate phonon frequency shifts it is important to 
compare the theoretical MVFF results with experimental data. Both the experimental 
acoustic and optical phonon shifts in Si nanocrystals are compared. 

Acoustic shift: The experimental values of acoustic phonon confinement (shown as dots 
with error bars in Fig. [3]) using Raman spectroscopy for Si nanocrystals embedded in a Si02 
matrix have been reported by Fujii et. aP^. The calculated acoustic phonon shifts in similar 
size free-standing circular [100] SiNWs using the MVFF phonon model compare quite well 
to the experimental data as shown in Fig. |3j providing frequency shifts between 20 and 30 
cm -1 . The MVFF calculated values show a power law of H^ -0 - 86 close to the results for [110] 
SiNWs (diamond dots in Fig. [3]) from RefP (W~ ' 95 ). The difference in the power laws may 
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arise from the wire surface passivation, structural relaxation and channel orientation. The 
theoretical predictions using the BCM modeP gives a cross-section dependence of W 
which overestimates the experimentally observed acoustic phonon shift (Fig. [3]). 

Optical shift: The optical phonon shift is experimentally measured using Raman spectroscopy' 1 
Figure [4] shows the experimental optical phonon shift result for Si nanorods (symbols) embed- 
ded in a Si02 matrix. With decreasing cross-section size the optical red-shift increases^. 
The theoretical optical phonon shift is calculated for free-standing [100] circular SiNWs 
using the MVFF model. The theoretical values of the optical shift agree well with the 
experimental data (Fig. [1]). The mismatch in phonon shifts at extremely small dimensions 
(W < 2nm) may be due to, (i) the structural relaxation in small cross-section wires^^ 
which is not taken into account in this study and (ii) the presence of S1O2 around SiNWs 
which might affect the results for smaller wires though the overlap of the optical phonon 
spectrum of these two materials is smalP^I. Hence, the MVFF model correctly captures the 
optical phonon shifts in SiNWs for W>2nm, however, for W < 2nm the MVFF phonon 
values without structural relaxation are not very accurate. Though the optical shift is very 
well captured by the Richter modeP^, the fitting parameters are chosen arbitrarily which 
removes the predictive capability from the model. 

In the next few sections we will discuss the structural properties of ultra-scaled SiNWs 
with different cross-section shapes and sizes. 

C. Acoustic Phonon shift in SiNWs 

Influence of cross-section shape and size: The acoustic phonon shifts are calculated using 
the MVFF model for different cross-section shaped SiNWs. All the wire geometries show 
a blue shift of the acoustic mode (this mode closely resembles the 'radial breathing mode' 
observed in perfectly circular wires^) as the wire cross-section size is reduced (Fig. [5] a). 
There is a clear cross-section shape dependence of the blue shifts where the blue shift follows 
the order Au lri < Au sq < Auh ex < Au C i r . The power exponent (d ac ) varies between 0.6 
and 0.856 (Table [ll|). The order of the pre-factor (A ac ) in Eq. §2§ value (Table |H| also shows 
the order of the acoustic blue shift for [100] SiNWs. 

With the increase in wire cross-section size two important things happen, (i) the geometric 
confinement on the phonons reduces, and (ii) the vibrational energy deviation between the 



surface and the inner atoms decreases (Fig. [5] First point explains the reduction in the blue 
shift with increase in cross-section size for all the SiNWs while the second point explains 
the order of blue-shift with shape. 

As the SiNW cross-section is reduced lots of low velocity phonon sub-bands start to 
appear above the acoustic branches (Fig. [2]) which are normally absent in bulk. These 
relatively flat bands cause phonon confinement. These zone center low energy phonon sub- 
bands cause acoustic blue shift. The surface atom vibrations become comparable to the 
inner core atom vibrations with size reduction. The blue shift is more pronounced when 
surface modes are less dominant. If the surface atoms vibrations are more than inner atoms 
then blue shift is suppressed^. 

The acoustic blue shift order can be explained by the difference between the vibrational 
energy density^ of the surface atoms and the inner atoms (AE). Surface atoms are the ones 
with dangling bonds. The spatial vibrational energy density can be calculated as^, 

E fT = J2 E [^K,) + O.5]^0 w -0^, (4) 

NA je[x,y,z] 

where NA is the total number of atoms in the unit cell, 'i' is the phonon sub-band, 'q' is 
the phonon wave- vector, 'j' is the polarization, Fbe is Bose-Einstein distribution for phonon 
population, is the phonon eigen frequency for sub-band T and wave-vector 'q' and 4>i,q,j 
is the phonon eigen vibration mode. 

Figure [5] b shows that the triangular wire exhibits the maximum difference between the 
surface and inner atom vibrational energy density for W = 2nm and 4nm, followed by the 
square, the hexagonal, and the circular wire. Figure [6] shows the spatial distribution of the 
energy density for the acoustic phonon branch for 2nm cross-section SiNWs. The surface 
atoms vibrate more intensively compared to the inner atoms. In circular and hexagonal 
nanowires all atoms vibrate with a similar intensity as indicated by the small AE value. 
However, in square and triangular wires the surface atoms vibrate more compared to the 
inner atoms leading to a larger AE value. 

The low frequency acoustic phonon modes are very similar to the breathing mode found 
in larger circular wires^. These eigen frequencies follow a 1/W d law (Fig. [5]), with d = 1 for 
a perfectly concentric circular wire^. The value of 'd' deviates from 1 as the cross-section 
shape deviates from the ideal circular geometry. 
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D. Optical Phonon shift in SiNWs 



Influence of cross-section shape and size: For optical phonon shift the MVFF model 



predicts a width exponent factor d^ between 1.8 and 2 (Table III) for all the cross-section 
shapes. The [100] triangular SiNWs show the maximum optical red shift (d opt = 1.795) 
whereas the square SiNWs show the minimum optical red shift (d op t = 2.06, close to the 
prediction in Ref.^) (Fig. [7]). Furthermore, with decreasing cross-section size (W) all the 
wires show an increasing optical red-shift (Fig. [7]). 



Many explanations have been reported in the literature for the optical red-shift in nanos- 
tructures like, the increasing surface-to-volume ratio (SVR)^^, compressive strain^, sur- 
face defects®! and optical phonon confinement) 21 * 23 * 2 -^) In the present study of free-standing 
SiNWs the observed shape and size dependence can be explained using the vibrational energy 
densitj^ G f the optical mode as well as the optical phonon group velocity (V^JP. 

Figure [8] shows the spatial distribution of the vibrational energy density for the optical 
mode in SiNWs with W = 4nm. The structure with the maximum spatial spread of the 
energy density has the most active optical modes^ since surface atom vibrations contribute 
to the optical modes. The symbols in Fig. [8] show the FWHM extent of the energy density. 
The square wire has the maximum energy spread while the triangular wire has the minimum 
energy spread which explains the highest red-shift in the triangular wires (least active optical 
modes). As the SiNW dimensions decrease, the energy spread also decreases which increases 
the red shift for all the shapes. 



Another reason for the red-shift can be attributed to the phonon group velocity. Figure [9] 
shows that the triangular wire has the smallest group velocity whereas the square wire has 
the highest group velocity. A smaller group velocity implies higher phonon confinement)^. 
Thus, the phonon group velocity gives the following order of optical phonon confinement; 
triangular > circular > hexagonal > square. Higher phonon confinement results in more 
red-shift. This again explains the observed shape dependence of optical red-shift in [100] 
SiNWs. 
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E. Orientation effect on phonon shifts 



Along with the shape and size of the SiNWs, the wire axis direction can also produce 
different phonon shifts. This indicates that the confined phonon dispersion is anisotropic. 
For brevity we present the orientation effects for square SiNWs only, however, these trends 
are also valid for other shapes. The acoustic phonon shift and the optical phonon shifts are 
shown in Fig. [10] a and b, respectively. The acoustic phonon blue shift exhibits a strong 



anisotropy (Fig. 10 a), however, the optical red shift has a very weak orientation dependence 



(10 b). This observed behavior of phonon shift can be explained by the phonon dispersion 
anisotropy in Si. In Si the maximum phonon frequency of the acoustic branch along [100], 
[110] and [111] directions are highly anisotropic ! 30 1 45 1 , whereas the maximum frequency of the 
optical branches are isotropic (~519.3 cm -1 ). The acoustic blue shift show the following 
behavior Aw^° > Aw^ 1 > Au^ . The cross-section size dependence for both kind of 



phonon shifts are provided in Table. IV The acoustic phonon shift exponent of [110] wires 
is 0.89 which compares reasonably to the exponent of 0.95 observed for [110] SiNWs in RefP 
(see Fig. 13b. 



F. Discussion 



The phonon spectrum shift in SiNWs with shape, size and orientation is useful in un- 
derstanding the experimental Raman shifts. The analytical expressions for the shifts can 
provide a good idea about the structure of the SiNWs which can be helpful in determining 
properties such as phonon-photon interactions. The evolution of the phonon shifts with 
increasing size also reflects the size (Wlim) at which the shape effects may be no more dis- 
tinguishable using Raman spectroscopy. From the expresssions of the phonon shifts (Eq. [2]) 
the Wlim value for acoustic blue shift comes out to be roughly 7nm and for optical shifts 
this value is lOnm. The analytical expressions are well suited for calculating the vibrational 
energy shifts depending on different physical aspects of the SiNWs. 



IV. CONCLUSION 



We have used a MVFF model to calculate the acoustic and optical phonon shifts in 
SiNWs. The theory agrees quite well with the experimental Raman shift data. At the 
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nanometer scale the acoustic phonon shift is quite sensitive to the (i) wire cross-section 
size and shape, and (ii) wire channel orientation whereas the optical phonon shift is mainly 
controlled by the wire cross-section shape and size. The different width exponents and 
the pre-factors from the acoustic and optical phonon shifts can become a useful tool in 
predicting the size and shape of the SiNWs. Triangular wires show the highest optical phonon 
confinement and the least structural symmetry. Analytical expressions for the shape and 
size dependence of the phonon shifts are useful for compact modeling of physical properties 
of ultra-scaled SiNWs. 
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TABLES 



TABLE I. Cross-section shape dependence of longitudinal and transverse sound velocity in W= 
3nm, [100] SiNW 



Shape 


V sn d,i(km/sec) 


V sn d,t{km/sec) 


Circular 


6.23 


4.50 


Hexagon 


6.18 


4.37 


Square 


6.26 


4.37 


Triangular 


5.79 


3.12 



TABLE II. Width parameters for acoustic phonon shift in [100] SiNWs 



Shape 


A ac (cm- 1 ) 


dac 


Circular 


140.05 


0.856 


Hexagon 


117.63 


0.824 


Square 


106.37 


0.81 


Triangular 


59.20 


0.6 



TABLE III. Width parameters for optical phonon shift in [100] SiNWs 



Shape 


A opt (cm,- 1 ) 


dopt 


Circular 


-70.49 


1.85 


Hexagon 


-62.59 


1.824 


Square 


-74.47 


2.06 


Triangular 


-106.05 


1.795 
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TABLE IV. Width parameters for acoustic and optical phonon shift in square SiNWs with different 
orientations 



Channel Or 


A ac {cm- 1 ) 


dac 


Aopt (cm- 1 ) 


dop 


[100] 


106 


0.81 


-74.62 


2.06 


[110] 


98.34 


0.89 


-66.47 


2.00 


[111] 


100.7 


0.87 


-57.17 


1.97 



FIGURE TITLES 




FIG. 1. Projected unit cell structures of free-standing [100] oriented silicon nanowires with (a) 
Circular, (b) Hexagonal, (c) Square, and (d) Triangular cross-section shapes. Width and height of 
the cross-section are defined using a single width variable W (width = height). These structures 
are at W = 2nm. 
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FIG. 2. Phonon dispersion in free-standing [100] oriented, W = 3nm, SiNWs with (a) Circular, (b) 
Hexagon, (c) Square, and (d) Triangular cross-section shapes. The first two branches (brown and 
black dot) are the flexural modest The next two branches (Acl, Ac2) are the acoustic branches^. 
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FIG. 3. Comparison of theoretical acoustic phonon frequency shift (line) obtained for circular 
[100] SiNWs using MVFF model with experimental phonon shift data (dots with error barsj^. 
The diamond dots show the abinitio based results for [110] SiNWs®. The BCM based results^ are 
shown using a dotted line. 
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FIG. 4. Comparison of theoretical optical phonon frequency shift (line), obtained for circular 
SiNWs using the MVFF model, and experimental data. The experimental optical phonon shift 
is for silicon nanorods embedded in SiO;P. The pre-factor and the exponent used for Richter's 
model is obtained arbitrarily by matching the experimental datsPl Though the match is good the 
fitting procedure lacks any predictive capability contrary to the MVFF model. 
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FIG. 5. (a) Effect of cross-section shape on the acoustic phonon shift in [100] SiNWs. (b) Difference 
in the average vibrational energy density of the inner and the surface atoms in 4 cross-section shape 
[100] SiNWs for W = 2nm and 4nm. 
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FIG. 6. Phonon energy density for the confined acoustic mode in 2nm x 2nm [100] SiNWs with 
(a) circular, (b) hexagonal, (c) square and (d) triangular shapes. Surface atoms vibrate more 
compared to the internal atoms. The difference in the energy of the surface and inner atoms is 
represented by AE. Triangular wires have largest AE which results in minimum acoustic phonon 
confinement (see Fig. [5]). 
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FIG. 7. Effect of cross-section shape on the optical phonon shift in [100] SiNWs. 
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FIG. 8. Spatial vibrational energy density for optical modes in different cross-section shaped 4nm 
x 4nm [100] SiNWs. The symbols show the full width half maximum (FWHM) of energy density 
which reflect how active the optical mode is. A larger spatial spread indicates more active optical 
modes. Optical modes in square wires are the closest to the bulk ones and triangular wires are the 
farthest. 




FIG. 9. Optical phonon group velocity in [100] SiNWs with 2nm x 2nm cross-section. 




FIG. 10. Effect of SiNW orientation on (a) acoustic phonon shift and (b) optical phonon shift. 
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